Illuminating ligand-induced dynamics in nuclear receptors through MD simulations.

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

How nuclear receptors transition between active and inactive forms: An energetic perspective.

Nuclear receptors regulate transcriptional programs in response to the binding of natural and synthetic ligands. These ligands modulate the receptor by inducing dynamic changes in the ligand binding domain that shift the C-terminal helix (H12) between active and inactive conformations. Despite decades of study, many questions persist regarding the nature of the inactive state and how ligands shift receptors between different states. Here, we use molecular dynamics (MD) simulations to investigate the timescale and energetic landscape of the conformational transition between inactive and active forms of progesterone receptor (PR) bound to a partial agonist. We observe that the microsecond timescale is insufficient to observe any transitions; only at millisecond timescales achieved via accelerated MD simulations do we find the inactive PR switches to the active state. Energetic analysis reveals that both active and inactive PR states represent energy minima separated by a barrier that can be traversed. In contrast, little or no transition is observed between active and inactive states when an agonist or antagonist is bound, confirming that ligand identity plays a key role in defining the energy landscape of nuclear receptor conformations.

Nuclear receptor interdomain communication is mediated by the hinge with ligand specificity.

Nuclear receptors are ligand-induced transcription factors that bind directly to target genes and regulate their expression. Ligand binding initiates conformational changes that propagate to other domains, allosterically regulating their activity. The nature of this interdomain communication in nuclear receptors is poorly understood, largely owing to the difficulty of experimentally characterizing full-length structures. We have applied computational modeling approaches to describe and study the structure of the full length farnesoid X receptor (FXR), approximated by the DNA binding domain (DBD) and ligand binding domain (LBD) connected by the flexible hinge region. Using extended molecular dynamics simulations (> 10 microseconds) and enhanced sampling simulations, we provide evidence that ligands selectively induce domain rearrangement, leading to interdomain contact. We use protein-protein interaction assays to provide experimental evidence of these interactions, identifying a critical role of the hinge in mediating interdomain contact. Our results illuminate previously unknown aspects of interdomain communication in FXR and provide a framework to enable characterization of other full length nuclear receptors.

Accelerating the discovery of alkyl halide-derived natural products using halide depletion.

Even in the genomic era, microbial natural product discovery workflows can be laborious and limited in their ability to target molecules with specific structural features. Here we leverage an understanding of biosynthesis to develop a workflow that targets the discovery of alkyl halide-derived natural products by depleting halide anions, a key biosynthetic substrate for enzymatic halogenation, from microbial growth media. By comparing the metabolomes of bacterial cultures grown in halide-replete and deficient media, we rapidly discovered the nostochlorosides, the products of an orphan halogenase-encoding gene cluster from Nostoc punctiforme ATCC 29133. We further found that these products, a family of unusual chlorinated glycolipids featuring the rare sugar gulose, are polymerized via an unprecedented enzymatic etherification reaction. Together, our results highlight the power of leveraging an understanding of biosynthetic logic to streamline natural product discovery.

The nuclear receptor LRH-1 discriminates between ligands using distinct allosteric signaling circuits.

Nuclear receptors (NRs) are transcription factors that regulate essential biological processes in response to cognate ligands. An important part of NR function involves ligand-induced conformational changes that recruit coregulator proteins to the activation function surface (AFS), ~15 Å away from the ligand-binding pocket. Ligands must communicate with the AFS to recruit appropriate coregulators and elicit different transcriptional outcomes, but this communication is poorly understood. These studies illuminate allosteric communication networks underlying activation of liver receptor homolog-1 (LRH-1), a NR that regulates development, metabolism, cancer progression, and intestinal inflammation. Using >100 µs of all-atom molecular dynamics simulations involving 74 LRH-1 complexes, we identify distinct signaling circuits used by active and inactive ligands for AFS communication. Inactive ligands communicate via strong, coordinated motions along paths through the receptor to the AFS. Activating ligands disrupt the "inactive" circuit and induce connectivity with a second allosteric site. Ligand-contacting residues in helix 7 help mediate the switch between circuits, suggesting new avenues for developing LRH-1-targeted therapeutics. We also elucidate aspects of coregulator signaling, showing that localized, destabilizing fluctuations are induced by inappropriate ligand-coregulator pairings. These studies have uncovered novel features of LRH-1 allostery, and the quantitative approach used to analyze many simulations provides a framework to study allosteric signaling in other receptors.

Two enhancer binding proteins activate σ54-dependent transcription of a quorum regulatory RNA in a bacterial symbiont.

To colonize a host, bacteria depend on an ensemble of signaling systems to convert information about the various environments encountered within the host into specific cellular activities. How these signaling systems coordinate transitions between cellular states in vivo remains poorly understood. To address this knowledge gap, we investigated how the bacterial symbiont Vibrio fischeri initially colonizes the light organ of the Hawaiian bobtail squid Euprymna scolopes. Previous work has shown that the small RNA Qrr1, which is a regulatory component of the quorum-sensing system in V. fischeri, promotes host colonization. Here, we report that transcriptional activation of Qrr1 is inhibited by the sensor kinase BinK, which suppresses cellular aggregation by V. fischeri prior to light organ entry. We show that Qrr1 expression depends on the alternative sigma factor σ54 and the transcription factors LuxO and SypG, which function similar to an OR logic gate, thereby ensuring Qrr1 is expressed during colonization. Finally, we provide evidence that this regulatory mechanism is widespread throughout the Vibrionaceae family. Together, our work reveals how coordination between the signaling pathways underlying aggregation and quorum-sensing promotes host colonization, which provides insight into how integration among signaling systems facilitates complex processes in bacteria.

The nuclear receptor LRH-1 discriminates between ligands using distinct allosteric signaling circuits.

Nuclear receptors (NRs) are transcription factors that regulate essential biological processes in response to cognate ligands. An important part of NR function involves ligand-induced conformational changes that recruit coregulator proteins to the activation function surface (AFS), ~15 Å away from the ligand binding pocket. Ligands must communicate with the AFS to recruit appropriate coregulators and elicit different transcriptional outcomes, but this communication is poorly understood. These studies illuminate allosteric communication networks underlying activation of liver receptor homolog-1 (LRH-1), a NR that regulates development, metabolism, cancer progression and intestinal inflammation. Using >100 microseconds of all-atom molecular dynamics simulations involving 69 LRH-1 complexes, we identify distinct signaling circuits used by active and inactive ligands for AFS communication. Inactive ligands communicate via strong, coordinated motions along paths through the receptor to the AFS. Activating ligands disrupt the "inactive" circuit by inducing connectivity elsewhere. Ligand-contacting residues in helix 7 help mediate the switch between circuits, suggesting new avenues for developing LRH-1-targeted therapeutics. We also elucidate aspects of coregulator signaling, showing that localized, destabilizing fluctuations are induced by inappropriate ligand-coregulator pairings. These studies have uncovered novel features of LRH-1 allostery, and the quantitative approach used to analyze many simulations provides a framework to study allosteric signaling in other receptors.

Identification of an Evolutionarily Conserved Allosteric Network in Steroid Receptors.

Allosteric pathways in proteins describe networks comprising amino acid residues which may facilitate the propagation of signals between distant sites. Through inter-residue interactions, dynamic and conformational changes can be transmitted from the site of perturbation to an allosteric site. While sophisticated computational methods have been developed to characterize such allosteric pathways linking specific sites on proteins, few attempts have been made to apply these approaches toward identifying new allosteric sites. Here, we use molecular dynamics simulations and suboptimal path analysis to discover new allosteric networks in steroid receptors with a focus on evolutionarily conserved pathways. Using modern receptors and a reconstructed ancestral receptor, we identify networks connecting several sites to the activation function surface 2 (AF-2), the site of coregulator recruitment. One of these networks is conserved across the entire family, connecting a predicted allosteric site located between helices 9 and 10 of the ligand-binding domain. We investigate the basis of this conserved network as well as the importance of this site, discovering that the site lies in a region of the ligand-binding domain characterized by conserved inter-residue contacts. This study suggests an evolutionarily importance of the helix 9-helix 10 site in steroid receptors and identifies an approach that may be applied to discover previously unknown allosteric sites in proteins.

Bile acids and the gut microbiota: metabolic interactions and impacts on disease.

Despite decades of bile acid research, diverse biological roles for bile acids have been discovered recently due to developments in understanding the human microbiota. As additional bacterial enzymes are characterized, and the tools used for identifying new bile acids become increasingly more sensitive, the repertoire of bile acids metabolized and/or synthesized by bacteria continues to grow. Additionally, bile acids impact microbiome community structure and function. In this Review, we highlight how the bile acid pool is manipulated by the gut microbiota, how it is dependent on the metabolic capacity of the bacterial community and how external factors, such as antibiotics and diet, shape bile acid composition. It is increasingly important to understand how bile acid signalling networks are affected in distinct organs where the bile acid composition differs, and how these networks impact infectious, metabolic and neoplastic diseases. These advances have enabled the development of therapeutics that target imbalances in microbiota-associated bile acid profiles.

Interactions governing transcriptional activity of nuclear receptors.

The key players in transcriptional regulation are transcription factors (TFs), proteins that bind specific DNA sequences. Several mechanisms exist to turn TFs 'on' and 'off', including ligand binding which induces conformational changes within TFs, subsequently influencing multiple inter- and intramolecular interactions to drive transcriptional responses. Nuclear receptors are a specific family of ligand-regulated TFs whose activity relies on interactions with DNA, coregulator proteins and other receptors. These multidomain proteins also undergo interdomain interactions on multiple levels, further modulating transcriptional outputs. Cooperation between these distinct interactions is critical for appropriate transcription and remains an intense area of investigation. In this review, we report and summarize recent findings that continue to advance our mechanistic understanding of how interactions between nuclear receptors and diverse partners influence transcription.

Ligand-induced shifts in conformational ensembles that describe transcriptional activation.

Nuclear receptors function as ligand-regulated transcription factors whose ability to regulate diverse physiological processes is closely linked with conformational changes induced upon ligand binding. Understanding how conformational populations of nuclear receptors are shifted by various ligands could illuminate strategies for the design of synthetic modulators to regulate specific transcriptional programs. Here, we investigate ligand-induced conformational changes using a reconstructed, ancestral nuclear receptor. By making substitutions at a key position, we engineer receptor variants with altered ligand specificities. We combine cellular and biophysical experiments to characterize transcriptional activity, as well as elucidate mechanisms underlying altered transcription in receptor variants. We then use atomistic molecular dynamics (MD) simulations with enhanced sampling to generate ensembles of wildtype and engineered receptors in combination with multiple ligands, followed by conformational analysis and correlation of MD-based predictions with functional ligand profiles. We determine that conformational ensembles accurately describe ligand responses based on observed population shifts. These studies provide a platform which will allow structural characterization of physiologically-relevant conformational ensembles, as well as provide the ability to design and predict transcriptional responses in novel ligands.

Differential Modulation of Nuclear Receptor LRH-1 through Targeting Buried and Surface Regions of the Binding Pocket.

Liver receptor homologue-1 (LRH-1) is a phospholipid-sensing nuclear receptor that has shown promise as a target for alleviating intestinal inflammation and metabolic dysregulation in the liver. LRH-1 contains a large ligand-binding pocket, but generating synthetic modulators has been challenging. We have had recent success generating potent and efficacious agonists through two distinct strategies. We targeted residues deep within the pocket to enhance compound binding and residues at the mouth of the pocket to mimic interactions made by phospholipids. Here, we unite these two designs into one molecule to synthesize the most potent LRH-1 agonist to date. Through a combination of global transcriptomic, biochemical, and structural studies, we show that selective modulation can be driven through contacting deep versus surface polar regions in the pocket. While deep pocket contacts convey high affinity, contacts with the pocket mouth dominate allostery and provide a phospholipid-like transcriptional response in cultured cells.

A phospholipid mimetic targeting LRH-1 ameliorates colitis.

Phospholipids are ligands for nuclear hormone receptors (NRs) that regulate transcriptional programs relevant to normal physiology and disease. Here, we demonstrate that mimicking phospholipid-NR interactions is a robust strategy to improve agonists of liver receptor homolog-1 (LRH-1), a therapeutic target for colitis. Conventional LRH-1 modulators only partially occupy the binding pocket, leaving vacant a region important for phospholipid binding and allostery. Therefore, we constructed a set of molecules with elements of natural phospholipids appended to a synthetic LRH-1 agonist. We show that the phospholipid-mimicking groups interact with the targeted residues in crystal structures and improve binding affinity, LRH-1 transcriptional activity, and conformational changes at a key allosteric site. The best phospholipid mimetic markedly improves colonic histopathology and disease-related weight loss in a murine T cell transfer model of colitis. This evidence of in vivo efficacy for an LRH-1 modulator in colitis represents a leap forward in agonist development.

Structural basis for an unprecedented enzymatic alkylation in cylindrocyclophane biosynthesis.

The cyanobacterial enzyme CylK assembles the cylindrocyclophane natural products by performing two unusual alkylation reactions, forming new carbon-carbon bonds between aromatic rings and secondary alkyl halide substrates. This transformation is unprecedented in biology, and the structure and mechanism of CylK are unknown. Here, we report X-ray crystal structures of CylK, revealing a distinctive fusion of a Ca2+-binding domain and a ß-propeller fold. We use a mutagenic screening approach to locate CylK's active site at its domain interface, identifying two residues, Arg105 and Tyr473, that are required for catalysis. Anomalous diffraction datasets collected with bound bromide ions, a product analog, suggest that these residues interact with the alkyl halide electrophile. Additional mutagenesis and molecular dynamics simulations implicate Asp440 in activating the nucleophilic aromatic ring. Bioinformatic analysis of CylK homologs from other cyanobacteria establishes that they conserve these key catalytic amino acids, but they are likely associated with divergent reactivity and altered secondary metabolism. By gaining a molecular understanding of this unusual biosynthetic transformation, this work fills a gap in our understanding of how alkyl halides are activated and used by enzymes as biosynthetic intermediates, informing enzyme engineering, catalyst design, and natural product discovery.

Structure and assembly of the diiron cofactor in the heme-oxygenase-like domain of the N-nitrosourea-producing enzyme SznF.

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates Nδ and Nω' of Nω-methyl-l-arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme's C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase-like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase-like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O2-reactive Fe2(II/II) complex and the structural basis for cofactor instability-a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon-hydrogen and carbon-carbon bonds to install olefins. Among ∼9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ∼25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

The Rho-GEF PIX-1 directs assembly or stability of lateral attachment structures between muscle cells.

PIX proteins are guanine nucleotide exchange factors (GEFs) that activate Rac and Cdc42, and are known to have numerous functions in various cell types. Here, we show that a PIX protein has an important function in muscle. From a genetic screen in C. elegans, we found that pix-1 is required for the assembly of integrin adhesion complexes (IACs) at borders between muscle cells, and is required for locomotion of the animal. A pix-1 null mutant has a reduced level of activated Rac in muscle. PIX-1 localizes to IACs at muscle cell boundaries, M-lines and dense bodies. Mutations in genes encoding proteins at known steps of the PIX signaling pathway show defects at muscle cell boundaries. A missense mutation in a highly conserved residue in the RacGEF domain results in normal levels of PIX-1 protein, but a reduced level of activated Rac in muscle, and abnormal IACs at muscle cell boundaries.

Allosteric regulation of thioesterase superfamily member 1 by lipid sensor domain binding fatty acids and lysophosphatidylcholine.

Nonshivering thermogenesis occurs in brown adipose tissue to generate heat in response to cold ambient temperatures. Thioesterase superfamily member 1 (Them1) is transcriptionally up-regulated in brown adipose tissue upon exposure to the cold and suppresses thermogenesis in order to conserve energy reserves. It hydrolyzes long-chain fatty acyl-CoAs that are derived from lipid droplets, preventing their use as fuel for thermogenesis. In addition to its enzymatic domains, Them1 contains a C-terminal StAR-related lipid transfer (START) domain with unknown ligand or function. By complementary biophysical approaches, we show that the START domain binds to long-chain fatty acids, products of Them1's enzymatic reaction, as well as lysophosphatidylcholine (LPC), lipids shown to activate thermogenesis in brown adipocytes. Certain fatty acids stabilize the START domain and allosterically enhance Them1 catalysis of acyl-CoA, whereas 18:1 LPC destabilizes and inhibits activity, which we verify in cell culture. Additionally, we demonstrate that the START domain functions to localize Them1 near lipid droplets. These findings define the role of the START domain as a lipid sensor that allosterically regulates Them1 activity and spatially localizes it in proximity to the lipid droplet.

Cutting in-line with iron: ribosomal function and non-oxidative RNA cleavage.

Divalent metal cations are essential to the structure and function of the ribosome. Previous characterizations of the ribosome performed under standard laboratory conditions have implicated Mg2+ as a primary mediator of ribosomal structure and function. Possible contributions of Fe2+ as a ribosomal cofactor have been largely overlooked, despite the ribosome's early evolution in a high Fe2+ environment, and the continued use of Fe2+ by obligate anaerobes inhabiting high Fe2+ niches. Here, we show that (i) Fe2+ cleaves RNA by in-line cleavage, a non-oxidative mechanism that has not previously been shown experimentally for this metal, (ii) the first-order in-line rate constant with respect to divalent cations is >200 times greater with Fe2+ than with Mg2+, (iii) functional ribosomes are associated with Fe2+ after purification from cells grown under low O2 and high Fe2+ and (iv) a small fraction of Fe2+ that is associated with the ribosome is not exchangeable with surrounding divalent cations, presumably because those ions are tightly coordinated by rRNA and deeply buried in the ribosome. In total, these results expand the ancient role of iron in biochemistry and highlight a possible new mechanism of iron toxicity.

Rewiring Ancient Residue Interaction Networks Drove the Evolution of Specificity in Steroid Receptors.

Understanding how changes in amino acid sequence alter protein dynamics and allosteric signaling would illuminate strategies for protein design. To gain insight into this process, we have combined molecular dynamics simulations with ancestral sequence reconstruction to explore conformational dynamics in two ancient steroid receptors (SRs) to determine how allosteric signaling pathways were altered over evolution to generate hormone specificity. In a broad panel of aromatized and non-aromatized hormones, we investigate inter-residue contacts that facilitate allosteric signaling. This work reveals interhelical interactions that act as ligand sensors and explain the structural and dynamical basis for ligand discrimination in SRs. These sensors are part of a conserved SR allosteric network and persist over long simulation time scales, indicating that evolutionary substitutions rewire ancient SR networks to achieve functional evolution. This powerful combination of computation, ancestral reconstruction, and biochemistry may illuminate allosteric mechanisms and functional evolution in other protein families.